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We introduce a new mean-field approach to the tt't" J model that incorporates both electron- like 
quasiparticle and spinon excitations as suggested by some experiments and numerical studies. It 
leads to a mean- field phase diagram which is consistent with that of hole and electron doped cuprates. 
Moreover, it provides a framework to describe the observed evolution of the electron spectral function 
from the undoped insulator to the overdoped Fermi metal for both hole and electron doping. The 
theory also provides a new non-BCS mechanism leading to superconductivity. 



The evolution of the electronic structure from the un- 
doped antiferromagnetic (AF) insulator to the overdoped 
metallic state of cuprates is a long standing problem. The 
plethora of anomalous behavior displayed by these mate- 
rials is particularly striking in hole underdoped samples, 
for which both experimental [1-6] and numerical [7-9] ev- 
idence suggests a dichotomy of the electronic excitations: 
excitations around the nodal points [k = (dz^,±^)] are 
well described as Landau's quasiparticles while those near 
the antinodal points [k — (7r, 0), (0, n)] show no signs of 
quasiparticle-like behavior. Some experimental [10] and 
numerical [7-9, 11] studies relate the absence of quasi- 
particles close to the antinodal points to the presence of 
excitations that only carry spin. 

In order to account for the aforementioned nodal- 
antinodal dichotomy, in this letter, a new mean- field 
(MF) approach to the tt't" J model is introduced which 
describes the low energy physics in terms of spinons and 
doped carriers. Spinons are electrically neutral fermions 
describing spin- 1/2 excitations. In the tt't" J model dou- 
ble occupancy is prohibited and the doped carriers cor- 
respond to the removal of a lattice spin, which inserts a 
unit charge and a spin- 1/2 in the system. Doped carriers 
are holes in the hole doped (HD) regime and electrons 
in the electron doped (ED) regime. For ease of speaking, 
below we refer to the doped carriers as dopons, which are 
spin- 1/2 charged fermions. We show that the new MF 
approach leads to a MF phase diagram that resembles 
the one of HD and ED cuprates. It also accounts for the 
doping evolution of the electronic structure, as seen by 
ARPES, in both HD and ED samples. 

We start with the 2D tt't" J Hamiltonian 
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where Uj = t,t',t" for first, second and third nearest 
neighbor (NN) sites and V projects out doubly occupied 
sites. The tJ model on-site Hilbert space, {| j) , \ [) , |0)}, 
includes states with either one or zero spin-| objects. 

To obtain the new MF theory we start with an enlarged 
on-site Hilbert space {|t0> , |j0> , |TT> , III) , lit) , III)} 
which contains either one or two spin-^ objects. The 
states 1 10), | jO) and the local singlet state -^=(|||) — | jf )) 



are the physical states that map onto the states ||), \[) 
and the vacancy state |0), respectively, in the tJ model 
on-site Hilbert space. The on-site triplet states, such 
as ^75 (IT!) + |!T))j are unphysical. We also introduce 
the fermionic representation for the first spin (the lat- 
tice spin), Si = ^fjcrfi, and the second spin (the doped 
spin), -(J-crdi, where a are the Pauli matrices. Here, the 

spinors f\ = [//^//j and d\ = [d\ +d\ .] are the spinon 
and the dopon creation operators on site i. Then, the 
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equals H t j in the physical Hilbert space and does not 
connect the physical and the unphysical sectors of the 
Hilbert space. H l enl is such that only local singlet states 
hop between different lattice sites whereas the unphysical 
local triplet states have no kinetic energy. Therefore, the 
dynamics included in a l enl effectively implements the lo- 
cal singlet constraint. The enlarged on-site Hilbert space 
contains at most one dopon. Hence, in (2), we intro- 



duce the projection operator V = Yli(^ 
which enforces the no double occupancy constraint for 
the d-fermion. By definition, the total number of dopons 
in the system equals the number of doped carriers. We 
are mostly interested in the low doping regime and, thus, 
below we drop the projection operators V in H^j 1 . 

The Hamiltonian H^j 1 is a sum of terms with up to 
six fermion operators. In the following we replace some 
multiple-fermion operators by their average so that the 
resulting MF Hamiltonian is quadratic in the operators 
f\ f, d) and and describes the hopping, pairing and 
mixing of spinons and dopons. 

The exchange Hamiltonian H J enl is decoupled by 
means of the d-wave ansatz [12] and becomes 
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bond and pairing MFs and ao is the Lagrange multiplier 
enforcing (fjfi) = 1. Also, J = (1 — x) 2 J, where we 
introduce the doping density x = (d\di). 

We now consider the hopping Hamiltonian H\ nl . Since 
the effective hopping amplitude of one hole in an AF 
background is renormalized by spin fluctuations, [13] 
we replace the bare £, t' and t" by the effective hop- 
ping parameters t\, ti and ts which are determined 
phenomenologically. The terms [(d\crdj).(iSi x Sj)] 
and (d\djSi.Sj) in H l enl are the sum of operators like 

A^j^ft-yfj^fj^hu and, in our decoupling scheme, 
only contribute to the MF spinon and dopon hopping 
terms. The first contribution comes from the averages of 
two d and two / operators ((d\ a dj^fJ 7 fj,s)) and yields 

the spinon NN hopping term l -f E(ij)ejvjv(/|/i + ft - c 0- 
The second contribution arises, instead, from taking the 
averages of the four / operators ({f} n fj,5f]^fi,v)), which 
reduce to (Si x Sj) and (Si.Sj), and adds up to the do- 
pon hopping term. We remark that, in the presence of lo- 
cal AF correlations, the vacancy in the quasiparticle state 
is surrounded by an AF-like configuration of spins. [9] To 
approximately account for this effect we assume that the 
spins encircling the vacancy in the one-dopon state are in 
a local Neel configuration. Therefore, we use (Si x Sj) = 
and (ASi.Sj) = (-1)^+^-^-^. Finally, to decouple 
the spinon-dopon interaction [(d\crdj).(Si-\-Sj)] we intro- 
duce b = (fjdi) and h = (^E^E^ nn fi d i+u), 
where u = ±£, u = ±x ± y and u = ±2x, ±2y for 
^ = 1,2,3 respectively. 

The resulting total MF Hamiltonian, written in terms 
of the Nambu operators rjj = faji^] = an d ^1 = 
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where at 



"(x^ _ ^x)(cosk x + cos/c y ) + a , ot x k = 
— ^A(cosk x — cosk y ), Pk — ^ [£i (cos fc^ + cosk y ) + 
2t 2 cos k x cos k y + ts(cos2k x + cos2/c y )] + &i and 7& = 
^2 cos fc^ cos fc^ + y (cos 2/^ + cos 2 fey) — /id, N is the lat- 

The 
and 



tice size and /i^ is the dopon chemical potential, 
eigenenergies of H^j F are ef k = ±VPfc 

±Vpk + V$k where 5 fc = /?£[(7fe + <) 2 + «) 
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ei 5 fc are the lowest, and 62, & the highest, energy bands. 

When 6o=^i=0 spinons and dopons do not mix. Then, 
the spinon sector of H^j F describes the same spin dy- 
namics as the slave-boson theory. [14] The dopon sector, 
on the other hand, determines the dynamics of doped 
quasiparticles. Here, the dopon only has intrasublattice 
hopping processes (see 7*.) due to the AF correlations in 
the spin average (Si.Sj) used to derive H^j F . In the HD 
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FIG. 1: (a) Regions in the (x-T) plane where A = or A / 
as well as bo = b± = or 60, b± 7^ 0. The dashed lines indicate 
the Trt described in the main text where long-range order in 
the dopon-spinon mixing channel is destroyed by vortex fluc- 
tuations, (b) The (x-T) phase diagram including the AF, SC, 
strange metal (SM), Fermi liquid (FL) and pseudogap with 
and without Nernst signal, labeled by N and PG respectively, 
regions. Both HD and ED cases are depicted in (a) and (b). 



regime, we choose £2 and ts so that 7^ approximately 
fits the high energy dispersion in ARPES data, [4, 15] 
which is isotropic around (f , f ) with a bandwith ~ 2 J 
for x w and whose gap at (71- , 0) closes around x ~ 0.3. 
[1] As a result, t% D = and tf? D = J -0.5 x f^J . In 
ED materials the electron pocket shows up around (7r, 0) 
instead [16] and we take t FD = 1.6 J 



0.6 x and 



t% u = 0.20 J + 0.3 x ^ J. In addition, choosing t x = 2 J 
correctly leads to the doping independent nodal disper- 
sion "kink" energy w ^ found in HD samples. [17] 

If spinons and dopons mix both bo and b\ are non- 
zero [18] and spin and charge dynamics become strongly 
coupled. Note that spinons and dopons are charge neu- 
tral and charged spin- 1/2 fermions while 60,1 are charged 
spin singlet fields. The condensation of 60,1 effectively at- 
tributes charge to spinons and, in the presence of spinon 
pairing (A 7^ 0), the system becomes superconducting 
(SC). Hence, the MF theory herein introduced provides 
a new route to the SC state via coherent spinon-dopon 
mixing or, equivalently, spinon-dopon pair condensation. 

The MF phase diagram in Fig. la contains four MF 
phases, all of which are observed in the cuprates: (a) 
d-wave SC when bo,bi,A 7^ 0; (b) Fermi liquid when 
bo,bi 7^ and A = 0; (c) pseudogap metal when 60, b\ = 
and A^0; (d) strange metal when bo,bi,A = 0. 

We note that the MF SC transition temperature is 
very high in the under doped regime. This is an artifact 
of the MF calculation since thermal fluctuations of the 
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phases of the condensates 60,1 are ignored. To crudely 
estimate the strength of 6o's phase fluctuations, we note 
that the NN electron hopping term in H t j induces a term 
— ^xZ^<ij>(^oi^oj + h.c). The resulting Kosterlitz- 
Thouless transition temperature Tkt = 0.9|ti|x&o> [19] 
above which the condensate average (bo) vanishes due to 
phase fluctuations, is plotted as the dashed-line in Fig. 
la. The state with long-range SC order only appears 
below Tkt (see Fig. lb). Above T^t, and in the un- 
der doped regime, there appear two distinct pseudogap 
metal regions marked by N and PG in Fig. lb. In re- 
gion N, which is located between the MF T c and Tkt, 
the non- vanishing magnitude of the MF order parameters 
60,1 leads to short-range SC correlations. This regime is 
observed experimentally, as suggested by the large Nernst 
signal measured in underdoped HD materials far above 
T c . [20, 21] In the PG region 6 ,i = and SC fluctuations 
become too small to be detected. 

In the above MF calculation we have ignored the 
AF phase. To include this state we further intro- 
duce the MF decoupling channels m = (— y x + % y (Sf ) 

and n = - (^=2,3 £«g«/ nn 4 a zdi + u + h.c.) 

that account for the staggered magnetization in the 
lattice spin and dopon systems respectively. We thus 
add 2J*Nm 2 - 47Vran - 2( J*m - n) Y, k ^l +(7r 7r) ^fc - 

2m Efc(7fc + / i rf) 7 ?l+( 7r , 7 r)^ to H u F \ where J* = A J and 
A = 0.31 is a renormalization factor that enforces the 
transition between AF and SC orders at x = 0.03 on the 
HD side. [22] Without addressing the issue of coexistence 
of AF and SC, we obtain the AF phase shown in Fig. lb. 
The hopping parameters in the ED regime favor intrasub- 
lattice hopping processes which do not frustrate AF. [23] 
Also, on the ED side dopons are located around (7r,0), 
which is away from the nodal points, thus weakening SC 
in the ED regime (it is destroyed at lower doping than in 
the HD regime). Therefore, AF order is very robust on 
the ED side where it extends over most of the SC dome 
and where it covers the pseudogap region N (in confor- 
mity with the lack of a vortex induced Nernst signal on 
these materials [24]). 

To compare the above MF theory to ARPES we note 
that c i<<7 = V^idl^jlj^f^fl^-dljl^h^V 
is the electron annihilation operator and the electron cre- 
ation operator in the HD and ED regimes respectively. 
Below, we ignore the incoherent contribution to the elec- 
tron spectral function and use c^ a = -^{d\ _ a + bofj_ a ) 
instead. Figs. 2a-2c show how the MF electron spectral 
function along the nodal direction evolves with hole dop- 
ing. These are self-consistent T = results concerning 
the SC phase. At T = only the two negative energy 
bands, namely e^ fc and e^~ fc , are occupied. For zero dop- 
ing the spectral function contains only a peak at e^ k 
(Fig. 2a). [25] Upon doping spectral weight is trans- 
fered from the e^ k to the k band, so that low energy 




FIG. 2: Electron spectral weights at T = 0. (a)-(c) Evolu- 
tion of the nodal direction electron spectral function with hole 
doping (top color scale). The white dashed line depicts the 
e ik band, (d)-(g) Electron spectral weight of the e^ k band 
states for different x in the HD regime (middle color scale). 
The white dashed line represents the minimum gap locus. 
The spectral weight at the node for x = 0.05, 0.12, 0.20, 0.25 
is 0.21, 0.35, 0.45 and 0.49 respectively, (h)-(i) Integrated 
electron spectral weight for x = 0.05,0.16 in the ED regime 
(bottom color scale). The energy window [— 0.15 J, 0.15J] 
was used. In (a)-(c) and (h)-(i) a Lorentzian broadening 
= -jq was used. 



quasiparticle weight develops above the parent insulator 
dispersion (hence inside the Mott gap!). As a result, in 
the underdoped regime two dispersive features arise. A 
linear dispersion crosses the Fermi level (lj=0) at a point 
that deviates from (§,§) toward (0,0). At higher en- 
ergy, a band that resembles the dispersion of the undoped 
AF samples carries most of the spectral weight. Remark- 
ably, this non-trivial behavior is also observed by ARPES 
[4, 15] and is to be contrasted with the conventional rigid 
band filling picture applicable to band insulators, namely 
that upon hole doping the chemical potential falls on top 
of the valence band forming hole pockets. 

Notably, Figs. 2d-2g show that spectral weight trans- 
fered from the high energy to the low energy band 
distributes in momentum space in agreement with the 
nodal-antinodal dichotomy displayed by ARPES data: 
(a) The spectral weight associated with each state in 
the k band develops on an arc-shaped region around 
the nodal direction. [4, 6] (b) In underdoped samples 
the spectral weight in the k band is depleted near the 
antinodal points and, as a result, the spectral structure 
in this fc-space region reflects only the high energy gap 
of (which is reminiscent of the AF insulator). [3] (c) 
The total spectral weight in the €[ h band increases with 
doping as the arcs extend to form a closed surface, (d) 
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The coherence peaks in the antinodal region only appear 
around and beyond optimal doping. [2] (e) A transition 
in the topology of the minimum gap locus from hole to 
electron-like at x « 0.20 is obtained. [5, 26, 27] 

Figs. 2h and 2i show that the MF low energy electron 
spectral weight distribution for the ED regime is also 
consistent with experiments. Indeed, at x = 0.05 there is 
AF order and an electron pocket is formed around (7r, 0) 
and (0,7r). [16] Further doping induces SC order and 
the d-wave SC quasiparticles develop spectral weight in 
the nodal region. As a result, a large "Fermi surface", 
ungapped only along the nodal direction, is observed in 
Fig. 2(i). [16, 28] 

To conclude, in this letter we introduce a new, fully 
fermionic, MF approximation to the tt't" J model. We 
also fit the MF parameters £1,2,3 to ARPES data to ar- 
gue that this MF approach is relevant to both HD and 
ED cuprates. As supported by the fact that £1,2,3 ~ J ^ 
the renormalization of the hopping parameters results 
from quantum spin fluctuations. Further work is required 
to properly understand the doping dependence of £1,2,3, 
which reflects the change of spin correlations as the sys- 
tem is doped. Remarkably, though, the MF approach 
correctly accounts for the evolution of low energy spec- 
tral weight from the undoped to the overdoped regime 
only by fitting the e^ k band to ARPES and by setting 
t\ = 2 J. [29] We stress that fitting the renormalized pa- 
rameters £1,2,3 to ARPES also leads to a relatively quan- 
titatively correct phase diagram. 

We analyze ARPES lineshapes in the cuprates in terms 
of a two-band description of the interplay between spin 
and charge dynamics. Related two-band interpretations 
were also proposed by numerical studies. [9, 30] In Ref. 
30 quantum Monte-Carlo results for the large U Hubbard 
model were interpreted in terms of two different states: 
(a) holes on top of an otherwise unperturbed spin back- 
ground and (b) holes dressed by spin excitations. Sim- 
ilarly, in Ref. 9 the nodal-antinodal dichotomy of the 
single hole tt't" J model was understood in terms of two 
types of states where: (a) the vacancy is surrounded by a 
staggered spin pattern and (b) the vacancy is surrounded 
by spins that screen the hole spin- 1/2 away. In this MF 
approach the doped carrier can also be surrounded by 
two different spin structures: (a) in the one-dopon state 
the vacancy is encircled by a local AF configuration of 
spins and (b) when the spinon and dopon mix the va- 
cancy is encircled by a local spin singlet configuration. 
[31] In case (a) NN hopping is strongly frustrated (in our 
MF approximation it is actually set to zero). However, in 
case (b) quasiparticles coherently hop between different 
sublattices - this fact shows up in the linear quasiparticle 
dispersion across the Fermi point [near (f , f )] (Figs. 2b- 
2c). The kinetic energy gain that follows the emergence 
of NN hopping stabilizes the formation of spinon-dopon 
pairs which lead to the SC phase. It also prevents the 
collapse of the chemical potential on top of the AF insu- 



lator band, [15] thus explaining the lack of hole pockets, 
in accordance with experiments. [4, 15] 
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